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ABSTRACT 

We describe a new approach for the determination of cosmological parameters 
using gravitational lensing systems with multiple arcs. We exploit the fact 
that a given cluster can produce multiple arcs from sources over a broad range 
in redshift. The coupling between the critical radius of a single arc and the 
projected mass density of the lensing cluster can be avoided by considering the 
relative positions of two or more arcs. Cosmological sensitivity appears through 
the angular size-redshift relation. We present a simple analytic argument for this 
approach using an axisymmetric, power-law cluster potential. In this case, the 
relative positions of the arcs can be shown to depend only upon the ratios of the 
angular-size distances of the source and between the lens and source. Provided 
that the astrometric precision approaches ~ 0.01 arcsec (e.g. via HST) and 
the redshifts of the arcs are known, we show that the system can, in principle, 
provide cosmological information through the angular-size redshift relation. 

We next consider simulated data constructed using a more general form 
for the potential, realistic sources, and an assumed cosmology. We present a 
method for simultaneously inverting the lens and extracting the cosmological 
parameters. The input data required are the image and measured redshifts 
for the arcs. The technique relies upon the conservation of surface brightness 
in gravitationally lensed systems. We find that for a simple lens model our 
approach can recover the cosmological parameters assumed in the construction 
of the simulated images. 
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1. Introduction 

One particularly promising approach for measuring the large-scale geometry of the 
universe is through the exploitation of gravitational lensing of distant galaxies and quasars. 
For some time gravitational lensing has been examined as a possible means of constraining 
cosmological parameters such as go and A. (e.g. [Turner, 1990|) To date most methods have 



focused on statistical arguments based on the cross sections for strong lensing of background 
objects (e.g. Cen et al., 1994 ). As an alternative we are investigating the feasibility of a 



method which makes more direct use of the dependence of gravitational lensing events on 
the geometry of the spacetime between the observer, the source and the lens. The method 
is motivated by recent images of cluster systems containing numerous arcs (e.g. Luppino 



et al., 1994] ; |Kneib et al., 1996| ). The method works by requiring consistency between both 



a single mass model for the cluster and a single angular-size redshift relation for all the 
lensed sources. 

We begin with an overview of gravitational lensing, followed by an analytic argument 
assuming an axisymmetric potential for the lensing cluster given by a power-law in density. 
We show that the astrometric precision afforded by HST is sufficient to distinguish between 
differing cosmologies through the angular-size redshift relation. We then investigate more 
realistic models for the lensing cluster and sources. However, in these cases the resulting 
arcs have more complicated morphology such that the lens and the geometry must be 
inverted via ray tracing and the minimization of a suitable merit function (e.g. one based 
upon the conservation of surface brightness and reasonable morphology of the sources, etc.) 
( |Kochanek et al., 1989 ). In section |3| we present an analysis of simulated data generated 



with an assumed form for the lensing potential and an arbitrary cosmology. With the 
simulated image and the redshifts of the arcs as input data we find that we can recover both 
the parameters describing the lensing mass distribution and the cosmology. In section |] 
we discuss the results of this analysis. We develop estimates of the bounds that could be 
set on cosmological parameters from a hypothetical system with the characteristics of our 
simulation, and we argue that this technique is a good complement to existing methods for 
the determination of cosmological parameters because it is purely geometrical rather than 
statistical in nature and is insensitive to evolutionary effects. In the concluding section we 
summarize the results of the previous sections, and we describe further refinements to our 
analysis which will be detailed in subsequent papers in this series. 
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2. Gravitational Lensing of Multiple Sources 

Gravitational lensing is sensitive to the gravitational potential of the lens, the angular 
position of the source object relative to the lens, and the angular size distances between 
observer, lens, and source (fig. |l]). This last is of principal interest, since the dependence 
of angular size distance on redshift has a direct dependence on cosmological parameters 
(e.g. |Sandage, 1988|) . Large arcs are produced by the near alignment of distant sources 



e. galaxies) and an intervening cluster of galaxies ( Paczynski, 1987 ; Soucail et al 



|1987t |Lynds fc Petrosian, 1989]) . For an axisymmetric (thin) lens with an on-axis source, 



the "critical radius" at which the source is observed depends only upon the projected 
surface mass density and the angular size distances to the lens and the source ([Blandford 



fe Narayan, 1992| ). Thus, by measuring redshifts for the gravitationally lensed sources 



and explicitly including the cosmological dependence of angular size distance on those 
redshifts, constraints could be set upon various cosmological models. We believe that 
the recently observed cluster systems with several multiply imaged arcs offer one of the 
best opportunities for obtaining meaningful constraints. We exploit the fact that sources 
distributed over a range in redshift can be imaged by the same lens. This should allow us 
to separate the effects of the lens mass distribution from the cosmological effects. 

The technique we propose involves inverting the lens system using methods similar 
to those presented in [Kochanek et al. (1989| ), but with one important difference; multiple 



sources at different redshifts are included. Inclusion of multiple source planes admits the 
possibility of cosmological dependence (q.v. §§ |2TT|J3.1| ), which is treated in our models 



by including the cosmological parameters as additional fitting parameters. For each 
cosmological model being tested the distances to all of the sources are computed using 
the redshifts and the particular cosmological model under consideration. The distances so 
calculated affect the quality of the fit, and this gives our method its cosmological sensitivity. 
While this does increase the number of free parameters in our model, we will show that the 
existence of multiple multiply imaged sources allows the modeled distances to be recovered 
with sufficient accuracy to provide interesting constraints on cosmological parameters. 



2.1. Analytic Treatment of an Axisymmetric System 

We can gain an understanding of the cosmological sensitivity of multiple-arc 
gravitational lens systems by considering a simple, if somewhat contrived, analytically 
treatable model. Accordingly, consider a cluster potential approximated as an axisymmetric 
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power law in density. 
The projected mass density is then 

Z M (R) = p a a b lR l -\ (2) 

where 

1= T d6cos h - 2 e, (3) 



which is independent of R. Integrating equation (fj) gives the total projected mass within 
the projected radius, R. 

M(R) = 27i Pa a b ^—R 3 - b . (4) 



For 6 = 2 this equation reduces to the familiar isothermal case (c.f. Preimer fc Sanders 



19921) . 



For simplicity, consider the case of on-axis sources, which would appear as "Einstein 
Rings" at the angular critical radii, 6e ■ Then 

2 _ AG N M(6 E ) D ls 

where Di, D s , and Di s are the angular size distances of the lens, source, and lens-to-source, 
respectively (c.f. |Schneider, Ehlers, fc Falco, 1992| , §2.1). Since R = D\Q we can combine 
equations (||) and such that: 

2 87lG N b 1 2 -b D lsp3-b f R \ 

0E ~ ~^P^ — h D l , (6) 

whence 

a -f 87TG " nn l> 1 D ^\^f Dls ^ (7) 

9e - {— paa —b Di ) IdtJ • (7) 

Thus, the angular radius of a given arc depends upon both the characteristics of the lens 
object and the angular-size distances to the source and lens-to-source. Note that for a 
particular lens system with multiple sources at differing distances the first expression on 
the left hand side of equation (0) is fixed, and only the second varies between sources. 
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Now consider the case of multiple arcs. The ratio of the Einstein radii for two sources, 
k and m is 
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We can combine this relation for three arcs, k, m, and n to eliminate b. 



In 



(fc) (m) 
D (k) D (m) 



In 



(9) 



The left side of equation (|9|) is an observational quantity which can in principle be measured 
directly, and the right side is a cosmological factor which can be calculated for any desired 
cosmological model, given the redshifts of the arcs. Thus, the relative radii of giant arcs has 
cosmological sensitivity and can provide a direct test of cosmological models. 

Standard error propagation techniques can be applied to equation (g) to determine the 
astrometric precision required to allow sensitivity to cosmological effects. Let the angular 
radii, 9 of the arcs be measurable with errors of ±er. Then the error in 



n = 



(10) 



is given by 
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We can evaluate equation ( [TT]) for a typical lens system, assuming a standard f2 = 1, A = 
cosmology. Consider a lens system (zi = 0.4) with three arcs (21,2,3 — 0.6,0.8,1.2). For 
9(z — > 00) = 30" the system produces arcs with angular radii #1,2,3 = V/8, 1 1'/8 , 15"7. These 
values in equation ( [LTD yield the result 



a-ji ~ O.lcr(arcsec). 



(12) 



The right hand side of equation @ can be evaluated for any desired cosmology. 
Figure shows calculated values of 1Z for two families of cosmological models: fiat 
cosmologies and A = cosmologies. From these we can conclude that a measurement error 
of g-ji ~ 0.001 would allow marginal discrimination between extreme cosmological cases. 
This corresponds to a measurement error in the positions of the arcs of cr arc ~ O'.'Ol. 
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We can compare the results of this section to the astrometric precision one might 
expect to obtain from presently available instruments. Examination of arcs imaged by HST 
shows that the cross section of a representative arc has a FWHM of (X'2. The arcs have 
a typical signal-to- noise ratio of 10. A simple monte carlo calculation indicates that the 
positions of the arcs should be measurable to an error of cr arc ~ 0'.'02. Such a measurement 
would result in marginal discrimination between cosmological models (fig. |2|). 

Real gravitational lens systems are expected to be more complex than these analytically 
soluble systems. Moreover, even for the simple potentials we have discussed, the off-axis 
sources present in real systems introduce an additional error term which can dominate the 
contribution to the total error. Consequently, the technique of comparing the astrometric 
positions of the observed arcs to their calculated positions for the models being tested is 
ill-suited for application to real systems. In section |3] we describe a numerical technique 
which addresses both of these issues. 



3. Evaluation of More Realistic Cases Using Simulated Data 

The surface brightness of a bundle of light rays is unchanged by gravitational deflection 
( Plandford fc Narayan, 1992p . We can exploit this property to construct a technique for 



analyzing systems more realistic than that described in section |2.1| . With this method 
comparing the surface brightnesses of image pixels mapped to the same point in the source 
plane replaces the measurement of arc positions as our primary test of the quality of a 
model. 



3.1. Lensing Potential 

The lensing potential is defined by 

a(0) = VV(0), (13) 

where ex is the reduced deflection angle depicted in figure [l], and is the position angle at 
which the source is observed. Thus, we see that the lens mapping is defined by 

(3 = 6- a 

(14) 

= 0-VV>(0). 



We adopt an elliptical lensing potential of the form (e.g. [Blandford fc Kochanek, 1987| ) 

2q 



V(x, y) = V + (1 + e c )x 2 + 2e s xy + (1 - e c )y 2 y. (15) 
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In this expression b is the angular critical radius (i.e. the radius of the critical curve), 
measured in arc-seconds, for an axisymmetric (e c = e s = 0) lens, s is the angular core 
radius, q is a function of the power-law index of the mass density distribution (0 < q < 0.5), 
and e c and e s are related to position angle, 0, of the major axis and the ellipticity, e, of the 
potential by: 

e c = ecos(20); (16) 
e s = esin(20). (17) 

It is worth noting that this form for the potential suffers from several defects. First, it 
retains the same ellipticity out to arbitrary distances; this is in contrast to the expected 
axisymmetric behavior at large distances from a mass distribution of finite extent. A 
related problem is that the surface mass density distributions produced by such a potential 
are unrealistic ("peanut-shaped") or unphysical (negative) at large radii. However, for 
sufficiently small e ( Blandford fc Kochanek ( 1987]) suggest e < 0.2) these anomalies occur at 



large enough radii that they do not affect the lensing behavior within the region of interest. 
While alternative potentials based on elliptical mass distributions have been studied, their 
lensing behavior is qualitatively the same within the region of interest as the elliptical 
potentials flKassiola fc Kovner, 19931 ). Consequently, we adopt the simpler potentials for 
this study. 

A potential of the form given in equation ( ]I5]) is sufficient to describe a system with a 
single source. However, for the case of multiple sources at differing redshifts lensed by the 
same mass distribution, b is different for each source. The constraint that each source is 
being lensed by the same mass distribution can be incorporated into a lensing potential of 
the form: 

1>(x, y) = Cfc^— V + (1 + e c )x 2 + 2e s xy + (1 - e c )y 2 )\ (18) 
2q 

where b is now the same for each source and is equivalent to the critical radius as z — > oo. 
For the fcth source (k is given by 

D {k) 

a = ^y- (19) 

With this form for the lensing potential the critical radius for each source is 

i 



b k = Cr q b. (20) 



Since the angular size distances D ls and are functions of f2 and A equation (|18D allows 
us to include the cosmological parameters as fitting parameters in our model, with the 
distances recalculated each time the cosmological portion of the model is varied. 
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3.2. Simulation of input data 

Our inversion requires two forms of input data: the redshifts of the lens and arcs 
and a high-resolution image of the system. In order to test the algorithm, simulated data 
sets were constructed using the following procedure. First a set of model parameters was 
selected. These parameters include the adjustable lens parameters in equation ([Tj|), the 
cosmo logical parameters, Q and A, the number of sources to be used, and the redshifts of 
the sources. The lens parameters were chosen to be representative of clusters of galaxies. 
The parameters used in our simulation are summarized in Table [l]. Source positions were 
arbitrarily selected but chosen such that strongly lensed arcs are produced. Realistic 
sources were constructed using a representative galaxy extracted from the Hubble Deep 
Field ( |Williams et al., 1996| , hereafter HDF). The redshifts for the lens and the four sources 



were chosen to be 0.2, 0.4, 0.6, 0.75, and 1.25, respectively. 

To incorporate HDF sources, several galaxies from the field were selected and tested; 
however, the models produced seemed to differ little, provided that the sources selected 
were scaled to a size appropriate for typical gravitational lensing sources. Ultimately, 
a single source was selected as being representative of the field and used for all of the 
sources in our simulations. Since the redshifts of most of the HDF sources are unknown, 
we scaled our source galaxy to produce simulated arcs of similar width to actual arcs. The 
background and peak intensities of these sources were also scaled in order to produce a 
desired signal-to-noise ratio in the arcs. 

Once model parameters are selected, each pixel in the image plane is mapped, via the 
lens mapping equation into each source plane using the angular-size redshift relation. 



Several methods were considered for selecting the appropriate surface brightness in the 
presence of several source planes, including the maximum brightness in any of the source 
pixels, the sum of the brightnesses of the source pixels, and the first (lowest redshift) source 
pixel encountered with a brightness above a particular threshold. In practice, since an 
image pixel does not generally map to bright regions in more than one source plane, there 
was little difference between the methods. Therefore, the surface brightness assigned to the 
image pixel is the maximum of the surface brightnesses of the source pixels to which it is 
mapped in each of the source planes. To produce a desired signal-to-noise ratio in the image 
a constant background level is added to each image pixel, and noise in the form of random 
Poisson deviates is added to the simulated image. The image produced by this simulation 
is shown in figure |^ (plate 00). 

The image produced by the simulation is stored for use as input to our inversion 
method. Also, the bright pixels in the image (i.e. the arcs) are identified with the source 
plane which produced them (i.e. the process fills the role of measuring the redshifts of real 
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arcs). These data, along with the redshifts of the source planes, form the input data for the 
inversion method. The model parameters used in the simulation are stored for comparison 
to the results of inversion attempts, but are not otherwise used. 

3.3. Merit function for inversion procedure 

A numerical model under consideration is evaluated by computing a merit function, /, 
for the model and comparing it to the merit functions of other models. Ideally, the merit 
function is so constructed that it has a unique minimum at the "true" solution. 

I ,1 

Kochanek et al. (1989) proposed a merit function based on the fact that surface 

brightness is conserved in gravitational lensing. We adopt a similar merit function, but 
with a few modifications which we find result in more robust behavior for our multi-source 
models. Our merit function is given by 

f = e + w + w s S. (21) 

In equation fl2"T| ) e, O and S are functions defined below which measure the properties that 
differentiate "good" models from "bad" . Specifically e measures the degree to which surface 
brightness is preserved in the lens mapping; O measures the degree to which the input 
image data is reproduced by that reconstructed from the model under consideration; and 
S measures the degree of fragmentation in the reconstructed sources (i.e. sources should 
resemble galaxies.) 

The quantities wo and wg are weights which ensure that the quantities being summed 
are all roughly of the same order of magnitude For the models presented here the weights 
used were 

w = 5, 
Ws = 0.005. 

The 'penalty' functions, O and S are designed to enforce physical constraints upon the 
numerical minimization algorithm. In practice the weight coefficients can be set to zero 
once the algorithm nears a solution. 

The components of the merit function are calculated as follows; an image of the 
system and measurements of the redshifts of the arcs visible in the system are required as 
input data. Arcs are identified as regions with a surface brightnesses above a particular 
"arc-threshold." The image plane is divided into regions corresponding to objects of 
measured redshift, blank regions, and regions containing blocking foreground or cluster 
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objects. At present we assume that redshifts can be measured for each arc in the system. 
This allows each arc to be identified with a particular source plane. Given a set of model 
parameters to be evaluated sources are reconstructed by ray tracing the image pixels back 
into the source planes corresponding to the appropriate redshifts. Each source pixel is 
then assigned the average of the values of the image pixels mapped to that source pixel. 
Additionally, multiply imaged pixels in each source plane are identified. Background pixels 
are not used in the reconstruction of the sources owing to the fact that they have no 
measured redshifts, and hence cannot be identified with a particular source plane. However, 
they are used in comparing the input and reconstructed images as described below. 

The function e measures the preservation of surface brightness and is defined by 

e = E _ n E W* - ( 22 ) 

i,j,k u ijV L ijk 1 ) m=l 

Here, k is summed over the source planes, and i and j are summed only over the multiply 
imaged pixels in each plane. is the value of the source pixel, is the multiplicity 
with which it is imaged, is the value of the mth image pixel identified with the source 
pixel, and is the noise in the pixel. 

In order to measure how well the image is reproduced, a reconstructed image is made 
from the source planes and the model parameters being evaluated. The function O is 
defined by 

= ^T,\ I a- R al (23) 

where refers to the surface brightness of pixel i,j of the input image, Rij refers to the 
surface brightness of the same pixel in the reconstructed image, and iV refers to the total 
number of pixels in the image. Since the sources are constructed to reproduce the observed 
arcs, the principal contribution to O is from pixels which were background in the input 
image, but which are brighter than the arc-threshold in the reconstruction. For real data, 
regions of the image for which there is no information (e.g. a blocking foreground object 
or non-transparent portion of the lens) would be excluded from the calculation of O. The 
sum is in practice limited to those pixels whose brightness exceeds the arc-threshold in 
the reconstructed image in order to prevent Poisson fluctuations in the background from 
dominating O. 

The morphology of galaxies at high redshift is not well known; however, we expect the 
surface brightness of neighboring pixels in a source plane to be correlated. Thus, we include 
a contribution to the merit function which measures the deviation of neighboring pixels in 
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the source plane. The function S is defined by 

S — — \Sijk S mn k\, (24) 

i,j,k m,n=i,j±l 

where the subscripted S refers to pixel values in the various source planes. Mathematically, 
S is the average difference between a pixel and its eight neighbors in the same source plane. 
This statistic measures the degree to which the source object is fragmented in the source 
plane; a contiguous source (e.g. one resembling a galaxy) is preferred over a fragmented 
source (e.g. one resembling a demagnified version of the arcs in the image). 

3.4. Determination of Model Parameters 

In order to extract the best fitting model parameters from the system, we must find the 
minimum of the merit function over nine parameters. Unfortunately, we were not able to 
design an automatic, globally convergent procedure for finding this minimum. The principal 
obstacle to creating such a procedure is that the merit function surface has much larger 
gradients in some principal directions than in others. Furthermore, far from the solution the 
merit function surface levels off to form a plateau with no large-scale gradient to indicate to 
an algorithm which way the minimum lies. The first problem can be overcome by making 
use of the underlying physical relationships among the parameters. The second problem can 
be addressed by noting that for each of the parameters there is an empirically determined 
"characteristic scale" within which an automatic minimization algorithm is convergent. 
Low-resolution scans can be made of the parameter space such that suitable initial starting 
values can be chosen. The characteristic scales for the model parameters were determined 
from an extensive search of the parameter space and are given in Table |l|. 

3.4-1- Computational Techniques 

To fully exploit the information contained within HST images of arc systems requires 
the use of a high-resolution grid of image and source pixels. The many evaluations of the 
merit function in a minimization algorithm places large demands on computer resources. 
The work presented here was performed using the SCAAMPQ collaboration's 64-processor 
Silicon Graphics Origin 2000. The facility has maximum throughput of 21 Gflops, although 
our calculations typically used only 16 of the 64 processors. The availability of this resource 



2 Scientific Applications on Arrays of Multiprocessors 



-12- 



allowed us to thoroughly study the parameter space in a reasonable amount of time. Details 
regarding the software developed for this procedure will be presented at a later date ( |Link 



A high-resolution scan of the parameter space was performed in order to determine 
the characteristic scales of the model parameters (Table [I]). We believe these scales will be 
generally applicable to gravitational lens systems parameterized in the form of equation (|I8D . 
Since the characteristic scales determine the size of the neighborhood of convergence, a scan 
of parameter space at this resolution can be used to select initial guesses for a minimization 
algorithm. Some representative slices of the parameter space are shown in figure |j. 

Having selected initial guesses we fix the cosmological parameters at reasonable values, 
the remaining parameters are determined using a standard downhill simplex technique 
( Press et al., 1988| ). Since the cosmological parameters have a small, second-order effect on 



the model we find that it is possible to arrive at a good solution for the lens parameters 
regardless of the initial guesses for cosmological parameters. 

Once the mass parameters for the system are obtained we solve for the cosmological 
parameters. Although small, the cosmological effects are strongly coupled only to the 
effects of the b parameter (c./. equation |18|.) Thus, with the exception of b, the lens mass 
parameters are held constant while the cosmological parameters are being varied. 

The correlation between b and the cosmological parameters can be understood as a 
tradeoff between the terms multiplying the front of the right side of equation ( jI5|) . This 
relationship can be used to assist the minimization algorithm in finding the cosmological 
solution once the mass parameters have been found. To further understand this relation, 
consider a simple system with only one source. Let us, however, describe this system using 
the form of the lensing potential given in equation ([181) . Then, all other things being equal, 
models for which 

£ll b 2(i-g) = cons t an t (25) 

will produce identical arc systems. In the case of multiple source models we find that the 
minimal surface can be approximated by 

£&2(i-«) = T ( 26 ) 

where T is a constant, and ( represents the distance ratio, Di s /D s for some sort of 
"characteristic" redshift of the sources. We find that by taking this characteristic redshift to 
be the median redshift of all of the sources, a reasonable approximation of the relationship 
between b and Q (with other parameters constant) can be obtained in the neighborhood of 
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the solution. However, since this approximation is strictly true only for the median redshift, 
we use equation ( p6[ ) only to provide an initial guess for a one-dimensional minimization in 
b at each evaluation of the merit function. 



4. Results 

We find that the minimization procedure outlined in the previous section converges to 
the global minimum of the merit function if the initial guesses for the model parameters 
are within a "neighborhood of convergence" about the correct values. This neighborhood is 
readily identifiable from the low-resolution parameter space scans described in section |3.4| . 
For a signal-to-noise ratio of 10 in the simulated data the lens parameters recovered by the 
algorithm are essentially indistinguishable from their fiducial values (q.v. Table p]), while 
the cosmological parameters recovered are accurate to within ~10%. 

We can gain some insight into the cosmological sensitivity of our technique by 
examining the morphology of sources and images reconstructed from a selection of 
cosmological models. In figures [5] and ^| (plates 00,00) we compare the simulated data 
with that reconstructed from both the best fit model and models with the cosmological 
parameters fixed at Q = 0, A = 1 and Q = 0.2, A = 0. In each of the latter two cases the 
non-cosmological parameters were varied in order to produce the best fit model for the 
cosmological model under consideration. Figure |5| shows a portion of the simulated image 
and corresponding regions of the image reconstructed from the model recovered by the 
algorithm and of images reconstructed from models with the cosmological parameters fixed 
at Q = 0, A = 1 and Q = 0.2, A = 0. Figure ^| shows the four source planes corresponding to 
each of these images. While the images and sources for the incorrect cosmological models 
are similar to the originals, they contain noticeable defects. In contrast, the reconstructions 
from the model arrived at by the minimization algorithm appear substantially better. 

Figure [j] shows contours of the merit function involving the cosmological parameters. 
Note that the cosmological parameters show strong correlations with b. This is due to 
the fact that these parameters affect the overall scale of the system through the (k in 
equation fllBD- Hence, for the first three panels b was adjusted to produce the minimal 
value of / at each point in the contour. This allows us to isolate the intrinsic effects of the 
cosmological parameters from those that result from changing the overall multiplicative 
constant in the lensing potential. 

Of particular interest is the slice in Q and A, since this characterizes our cosmological 
sensitivity. For our particular simulation the cosmological sensitivity is primarily to fl 
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rather than A. This is partly a result of the fiducial model we have chosen for this study 
and partly a consequence of the range of redshifts chosen for the simulation. By choosing 
Q = 1, A = 0, a strongly mass-dominated model, for our fiducial cosmology, we effectively 
mask the influence of A unless A is comparable to Q. Moreover, the angular size-redshift 
relation, upon which this technique depends for its cosmological sensitivity, is sampled only 
over the range redshifts of spanned by lens and the sources. Consequently, cosmological 
models which produce similar angular size-redshift relations over those redshifts will "look" 
similar. Examination of equation fll8|) reveals that the cosmological parameters enter into 
the lensing potential through the leading multiplicative term: 

t = C(*) (27) 

Figure || illustrates the dependence of this expression on redshift for selected cosmological 
models. For each cosmological model we evaluate the quantity £ in equation (p7|) , using 
the values of b and q which minimize /, and compare this to the same quantity evaluated 
for the fiducial cosmology. The differences are plotted as a function of redshift for each 
cosmological model. For comparison we show three cosmological models of interest: model 
I, a flat, low-omega model; model II, a A = 0, low-omega model; and model III, a model 
along the near-minimal surface of / in the Q-X plane. Note that the residuals for model 
III are substantially smaller than those for the other models over the range of redshifts 
represented in the simulated data. Thus, it is understandable that model III is a better fit 
(as measured by its relatively low value of /) than models I and II. Therefore, the shape 
of the merit function contours in the O-A plane results from the behavior of equation (27) 
over the range of redshifts sampled by the lens and sources. Specifically, the contours are 
elongated in the direction of cosmological models that produce values of £ similar to those 
of the fiducial model at the sampled redshifts. 



4.1. Monte carlo estimate of confidence regions 

We have identified four principal sources of error in the parameters determined using 
our method. These are: the noise in the input image, errors in the redshift measurements, 
deviation of the actual lens potential from our parameterized form, and any density 
fluctuations along the line of sight which alter the angular size-redshift relation. We 
consider only the first two sources of error in this study; although the others will be 
addressed in subsequent papers in this series. 

In order to evaluate these first two sources of error, we have performed monte carlo 
calculations to estimate confidence regions in the cosmological parameters. At each iteration 
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of these calculations new simulated data are constructed from the fiducial model described 
in the previous sections, using a new seed to create an independent image of the system 
with the same signal-to-noise ratio as the original. Next, errors with a Gaussian distribution 
are added to the redshifts in the model, and the merit function, /, is calculated for the 
fiducial parameters and the perturbed redshifts. The process is repeated until a sufficiently 
large sample is collected; in this case 100 samples were used. The contour in the Q-X 
parameter space corresponding to the 68.5 percentile of / is adopted as an estimate for 
the la confidence region for Q and A. Figure |] shows the la confidence regions obtained 
through this technique for several values of a z , the standard error in the redshifts. With a 
peak signal-to-noise ratio of ten in the image arcs, the contribution from the measurement 
errors in the redshifts is by far the dominant effect. 

If we adopt a z = 0.001 as a typical error in current measurements the allowed region 
of the Q-X parameter extends over a broad range in A and a modest range in Q. However, 
much of the parameter space is nevertheless excluded. In particular, our method strongly 
strongly discriminates between "interesting" cosmo logical models, viz. flat cosmologies 
(Q + A = 1) and A = cosmologies. For these classes of models the constraints from this 
approach are strong. Results for these specific cases are summarized in Table |[ It is worth 
noting that the allowed range of cosmological parameters is significantly reduced if the error 
in the arc redshifts can be reduced to a z w 0.005. While such measurements are difficult, 
they are not beyond the reach of 10-m class telescopes. 



5. Concluding Remarks 

Our inversion technique is able to extract both the lens and the cosmological parameters 
used to construct the simulated data. Moreover, the confidence regions obtained in our 
monte carlo calculation suggest that physically interesting constraints could be set on 
cosmological parameters through this method. The requisite data, high-resolution images 
from HST and redshifts for the arcs, are are currently obtainable. We are encouraged by 
recent observations which include systems with multiple sources QLuppino et al., 1994j ; 
Kneib et al., 1996| ), and by the fact that some systems appear to have sources at several 



redshifts ( [Hogg et al., 1996] ). 



We have made several simplifying assumptions in our analysis, particularly the use of 
the same analytic potential for both simulating and analyzing the data. We are currently 
engaged in further work studying the cosmological constraints which result from applying 
this technique to a potential whose form is not known a priori, and which may contain 
substructure. Presumably these factors will introduce additional error which may weaken 
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the constraints which can be set on cosmological parameters using this technique. However, 
in view of the success we have experienced with our simulations to date we believe that 
the errors so introduced will be surmountable and that this method is capable of placing 
interesting bounds on cosmological parameters. 
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Fig. 1. — The geometry of gravitational lensing in the thin lens approximation. For a light 
ray passing through the lens plane at impact parameter £, the reduced deflection angle (a) is 
given by the gradient of the lensing potential. The undeflected position of the source (f3) is 
related to the observed position (9) by the lens mapping (equation [14] in text). All distances 
are angular size distances. 
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Fig. 2. — The dependence of relative angular radii of gravitational arcs on cosmological 
parameters for an analytically soluble axisymmetric lens model (see § |2.1|) . Here 1Z = 
In f |H / In (J±j . Computed values of 7Z are shown for two families of cosmological models. 
For comparison the we include a conservative estimate of the uncertainty for measurements 
obtainable using HST. 
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Fig. 3. — The simulated image used in this study. The inset shows a source extracted from 
the Hubble Deep Field and used to construct the simulated data. The source was scaled and 
resampled to produce sources at each redshift in the simulated system. 
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Fig. 4. — Contours of the merit function, /, for a selection of 2-D slices in the parameter 
space. In each case the cross marks the fiducial parameters used in the simulation. Note the 
correlation between the power-law index of the lens potential, q, and the overall scale factor, 
b (q.v. equation p~8|). 
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Fig. 5. — A comparison of a fixed region within (clockwise from upper-left) the simulated 
image, the image reconstructed from the best fit model (see also Table H), and the images 
reconstructed from the best fit models with cosmological parameters fixed at O = 0, A = 1, 
and Q = 0.2, A = 0, respectively. Note that the best fit model reproduces the simulated 
image well, while the others have noticeable defects. 
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Fig. 6. — A comparison of the input data and reconstructed source planes for the models 
shown in figure [5]. From top to bottom are the four source planes corresponding to the 
redshifts of the arcs in the system. From left to right are the input data, the best fit model, 
and the best fit models with Q = 0, A = 1 and Q = 0.2, A = 0. Note that as was the case for 
the reconstructed images, the best fit model reproduces the simulated data well, while there 
are noticeable defects in the sources corresponding to the other models, particularly at the 
two redshift extremes. 
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Fig. 7. — Contours of the merit function, /, for parameter space slices including the 
cosmological parameters (f2 and A). In slices not involving b, we adopt the value of b which 
minimizes / with other parameters fixed. In each case the cross marks the parameters used 
in the simulation. Note in panel (d) the correlation between b and Q which results from the 
fact that both parameters affect the overall scale of the system. Note also the existence of a 
well-defined minimum in the cosmological parameters. This demonstrates the ability of this 
technique to discriminate between cosmological models. 
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Fig. 8. — Residuals of the quantity £ in equation ^ with respect to the fiducial model as 
a function of redshift. The crosses mark the redshifts of the arcs in our simulated data. 
Models which produce small residuals at these redshifts should produce "better" fits with 
this technique. Model I is a flat, low-omega model, model II is a A = 0, low-omega model, 
and Model III is a model along the curve of the elongated contours of / in the Q-\ plane 
(figures [7| and |). Note that model III, which produces the best fit of these three, also 
produces the smallest residuals of £ with respect to the fiducial model. 
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Fig. 9. — A contour of the merit function, /, in the Q-X subspace showing the confidence 
region estimated for several values of a z using the procedure in section 4.1 . For o z = 
(i.e. no measurement errors in the redshifts of the arcs) the estimated confidence region is 
approximately the size of the cross which marks the fiducial model. Although the allowed 
region is fairly large with typical redshift errors, along the directions of astronomical interest 
(Q + A = 1, shown by the dashed line, and A = 0) the constraints are quite strong (q.v. 
Table 0). 
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Table 1. Simulation Parameters 



Parameter Characteristic Length Fiducial Value' 





large b 


1.0 


X 


large b 


0.0 


b 


r/5 


15"0 


Q 


0.05 


0.45 


s 


1" 


2'/5 


e c 


0.05 


0.025 


e s 


0.05 


0.025 


x 


1" 


070 


Vo 


1" 


O'/O 



a The values of the parameters recovered by this method (see text) 
are indistinguishable from their fiducial values. 

b "Large" means that the characteristic length of the parameter 
extends to the entire physical range of the parameter. 
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Table 2. Cosmological Bounds vs. Redshift Errors. 



Case 




a z = 0.001 


a z = 0.0005 


n + x = 


1 


n > 0.8, A < 0.2 


n > 0.9, A < 0.1 


A = 





0.8 < n < 1.2 


o.9 < n < i.i 


n = 


i 


A < 0.8 


A < 0.5 



